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ABSTRACT 

We present high-contrast observations of the circumstellar environment of the Herbig Ae/Be star 
HD100546. The final 3.8 fim image reveals an emission source at a projected separation of 0.48"±0.04" 
(corresponding to ^47±4 AU) at a position angle of 8.9°±0.9°. The emission appears slightly extended 
with a point source component with an apparent magnitude of 13.2 ± 0.4 mag. The position of the 
source coincides with a local deficit in polarization fraction in near-infrared polarimetric imaging data, 
which probes the surface of the well-studied circumstellar disk of HD 100546. This suggests a possible 
physical link between the emission source and the disk. Assuming a disk inclination of ^47° the 
de-projected separation of the object is ^68 AU. Assessing the likelihood of various scenarios we favor 
an interpretation of the available high-contrast data with a planet in the process of forming. Follow- 
up observations in the coming years can easily distinguish between the different possible scenarios 
empirically. If confirmed, HD100546 "b" would be a unique laboratory to study the formation process 
of a new planetary system, with one giant planet currently forming in the disk and a second planet 
possibly orbiting in the disk gap at smaller separations. 

Subject headings: stars: formation — planets and satellites: formation — protoplanetary disks - 
planet-disk interactions — stars: individual (HD100546) 



1. INTRODUCTION 

To possibly extend the ongoing census of exoplanet 
demographics from our solar neighborhood to the en- 
tire Milky Way we need to understand planet formation 
and its dependence on initial physical and chemical con- 
ditions. A key step is to directly detect and charac- 
terize forming planets in their natal environment. Re- 
cently, based on sparse aperture masking observations, a 
few low-mass companion candidates have been revealed 
in the gap of their host star's transitional disks (e.g., 
LkCal5 b, T Cha b; Kraus & Ireland 2012; Huelamo 
et al. 2011), but in some cases scattered light from the 
disk rim or other disk structures may be a valid expla- 
nation for the observed features (e.g., Cieza et al. 2013). 
Until now no protoplanet has yet been found embedded 
in the optically thick gas and dust disk of its host star. 
Here we present high contrast imaging data revealing a 
protoplanet candidate embedded in the disk around the 
Herbig Ae/Be star HD100546. 

HD100546 (see Table 1 for stellar properties) has a 
complex circumstellar environment consisting of an inner 
disk from ^0.2-4 AU, a disk gap from ~4-13 AU, and a 
large outer disk from ~13 AU out to a few hundred AU 
(e.g., Benisty et al. 2010; Tatulli et al. 2011). A massive 
planet was suggested to be orbiting in the gap (Bouwman 
et al. 2003). Based on asymmetries in the line profile of 
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TABLE 1 
Basic parameters of HD100546. 



Parameter 


HD100546 


Reference a 


RA (J2000) 


ll h 33 m 25 s .44 


(1) 


DEC (J2000) 


-70°11'41".24 


(1) 


J 


6.43±0.02 mag 


(2) 


H 


5.96±0.03 mag 


(2) 


K s 


5.42±0.02 mag 


(2) 


Mass 


2.4±0.1 M 


(3) 


Age 


5...M0 Myr 


(3),(4) 


Distance 


97t 4 P c 


(1) 


Sp. Type 


B9Vne 


(5) 



a (l) van Leeuwen (2007), (2) 2MASS point 
source catalog (Cutri et al. 2003), (3) van 
den Ancker et al. (1997), (4) Guimaraes et al. 
(2006), (5) Houk & Cowley (1975). 

[OI] (Acke & van den Ancker 2006) and OH (Liskowsky 
et al. 2012) dynamic evidence for such an object was 
found. In the case of OH the emission is thought to arise 
from an eccentric inner rim of the outer disk with the ec- 
centricity being introduced by a planet. The outer disk 
has been resolved at multiple wavelengths including scat- 
tered light (e.g., Augereau et al. 2001; Grady et al. 2001; 
Ardila et al. 2007) where it shows peculiar features such 
as large-scale spiral arms. The remaining disk mass is es- 
timated to be 10" 2 - 10" 3 Mq (e.g., Panic et al. 2010). 
Recently, polarimetric differential imaging (PDI) in the 
near-infrared (NIR) revealed distinct sub-structures in 
the innermost few tens of AU of the disk (Quanz et al. 
2011b). 

2. OBSERVATIONS AND DATA REDUCTION 

HD100546 was observed with VLT/NACO (Lenzen 
et al. 2003; Rousset et al. 2003) and its Apodozing Phase 
Plate (APP) coronagraph (Kenworthy et al. 2010), which 
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TABLE 2 

Summary of deep imaging observations in pupil tracking mode. 



Parameter 


TJTV1 nn^ylfi 


riJJ_LUU04o 




I— I f>rmonnof^ l 


1— T otvi i annovo O 


Date 


2011-05-30 


2011-07-13 


UT start 


22h:48m:18s 


23h:00m:49s 


UT end 


00h:14m:58s 


00h:26m:lls 


NDIT x DIT a 


200 X 0.15 s 


200 X 0.15 s 


NINT b 


130 


126 


Parallactic angle start 


-17.33° 


40.36° 


Parallactic angle end 


10.33° 


63.93° 


Airmass range 


1.45. . . 1.43 


1.54. . . 1.75 


Mean DIMM seeing [A=500 nm] 


0.6" 


0.9" 


(To)mean / (To)min / (lj)max C 


4.4 / 0.0 / 9.7 ms 


1.9 / 0.0 / 4.0 ms 


PA start d 


-108.35° 


129.34° 


PA cnd d 


-80.86° 


153.04 ° 



a NDIT = Number of detector integration times (i.e., number of individual 
frames); DIT = Detector integration time (i.e., single frame exposure time). 
b NINT = Number of data cubes. 

c Average, minimum and maximum value of the coherence time of the 
atmosphere in data cube. Calculated by the Real Time Computer of the 
AO system. 

d Position angle of camera adaptor. 



was already used in earlier exoplanet imaging projects 
(Quanz et al. 2010, 2011a; Kenworthy et al. 2013). We 
used the L27 camera (~ 27 mas pixel -1 ) with the L' 
filter (A c = 3.8 /im, AA = 0.62 /im) in angular differ- 
ential imaging mode (ADI; Marois et al. 2006). Two 
datasets were taken on two different nights, referred to 
as "hemisphere 1" and "hemisphere 2", with a 180° offset 
in position angle of the APP to allow its high-contrast 
half to cover most of the circumstellar environment. Us- 
ing the "cube mode" option, all exposures, each 0.15 s 
long, were saved individually. The core of the stellar 
PSF was slightly saturated. Data cubes consisting of 
200 exposures were taken using a 3-point dither pattern 
along the detectors x-axis with roughly 7" separation. 
Table 2 summarizes the observations and the observing 
conditions. 

For photometric calibration we observed HD100546 in 
the NB3.74 filter (A c = 3.74 /im, AA = 0.02 /im; cen- 
tered on the Pf 7 line) . The observing strategy was iden- 
tical to the deep science observations, but these expo- 
sures were in the linear detector regime (0.2 s exposure 
time) and only 2 data cubes with 150 exposures each 
were recorded. No photometric or astrometric standard 
star was observed. 

The basic data reduction steps (bad pixel cleaning, sky 
subtraction, image alignment) were done in a similar way 
as described in Quanz et al. (2010). During the alignment 
process the images were re-binned to twice their original 
resolution. From each image in the stack of aligned ex- 
posures we created a 2"x2" sub- image centered on the 
star. Individual images showing bad AO correction or 
not covering the full size of the sub-images were disre- 
garded. In the end we had a stack of 16,117 images for 
hemisphere 1 and 18,916 images for hemisphere 2. 

The PSF-subtraction was done using the principle 
component analysis based software package PynPoint 
(Amara & Quanz 2012), and the results were confirmed 
by using the LOCI algorithm (Lafreniere et al. 2007). 
For the final PynPoint images we used 80 PCA coeffi- 
cients and kept only the best 12000 images in terms of 
total residuals over the whole image frame. After PSF 



subtraction, each image was de-rotated to the same field- 
orientation and we computed the mean image of the im- 
age stack clipping data points that were beyond 2.5cr of 
the mean. The results shown below are robust against 
all of these numbers (see next section). 

For LOCI we median combined 20 consecutive expo- 
sures into a single image in those data cubes where this 
was possible, resulting in 734 stacked images. The LOCI 
algorithm was then applied to this stack of images using 
the following LOCI parameters: FWHM=8 px, ^=0.75, 
dr=8, Na— 500. The choice of these values reflect the fact 
that the images have been re-binned to twice their origi- 
nal resolution (see above). All final images (see, Figure 1) 
were smoothed with a circular gaussian with a width of 
3 pixels. 

3. RESULTS AND ANALYSIS 

3.1. Detection of an emission source 

In Figure 1 an emission source is revealed north of the 
central star in the hemisphere 1 dataset. To examine the 
robustness of this detection we did a series of tests: 

• We varied the number of PCA coefficients used in 
PynPoint (between 20 and 200). 

• We split the dataset in half in different ways in- 
cluding random selections of ^50% of the images. 

• We applied the LOCI algorithm as independent re- 
duction approach. 

In all cases we found a bright feature at the same loca- 
tion. In addition, we analyzed the hemisphere 2 dataset 
in exactly the same way without finding a persistent 
source at any location in the final image. 

Using the approach described in Quanz et al. (2011a) 
and a 8-pixels wide aperture (~ 1-FWHM), the source 
has a signal-to-noisc of ~15 in the final LOCI image. 
The detected emission appears slightly elongated in the 
northern direction (Figure 1). To estimate the photome- 
try and astrometry of the source we did a detailed anal- 
ysis inserting fake negative planets (details see below). 
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Fig. 1. — NACO/APP L' images of the circumstellar environment of HD100546. From left to right: Final PynPoint images of hemisphere 
1 and hemisphere 2 and final LOCI image of hemisphere 1. An emission source is clearly detected in left and right panel. The shaded area 
indicates the region that was only covered by the low sensitivity hemisphere of the APP. The images are scaled with respect to their peak 
flux. 



It showed that the observed emission can be explained 
with a point source plus some extended component. For 
the point source the projected separation amounts to 
~0.48"±0.04" (~ 47±4 AU). The uncertainties in the ex- 
act location of the central star and the point source com- 
ponent in x and y on the detector were 0.5 and 1 pixel, 
respectively. The position angle of the point source with 
respect to the star is ~8.9°±0.9°. This error excludes 
any systemic error in the orientation of the camera with 
respect to the true celestial north, which is estimated to 
be <0.5° based on calibration data from an ongoing large 
imaging program (PI: J.-L. Bcuzit). 

To estimate the brightness of the point source we in- 
serted artificial negative planets in the individual expo- 
sures and re-ran PynPoint. For the fake objects we 
used an unsaturated PSF of HD100546 from one of the 
photometric calibration datasets. To scale the flux of 
these objects, the difference in exposure time between 
the science and the calibration images had to be con- 
sidered as well as the transmission curves of the two 
different filters 6 . Using published L-band spectra for 
HD100546 from ISO and VLT/ISAAC (Geers et al. 2007, 
and references therein) we derived a throughput fraction 
of ^0.074 ± 0.002 for the narrow band filter compared to 
the broadband L' filter. The error arises from changes 
in the Pf 7 line emission in HDI00546 between the two 
datasets suggesting that the NB3.74 filter traces variable 
accretion activity. Also the whole NIR and MIR contin- 
uum varies with an offset of a factor of ~1.25 between 
the two datasets of Geers et al. (2007), which impacts 
the error in the final photometry. 

For our contrast estimates we varied the brightness of 
the injected fake sources in steps of 0.1 mag and used 
two independent methods. First, we searched for a neg- 
ative source that, when subtracted, yielded a remaining 
flux at the object's location similar to the flux level in 
the surroundings, i.e., in the extended flux component. 
Secondly, we canceled out all the flux at the object's lo- 
cation. Our best estimate contrasts were AL' = 9.0 ±0.3 
mag and AL' = 8.3 ± 0.3 mag for both methods, respec- 
tively. A contrast closer AL' w 9.0 mag appears to be 
more likely, as strong residuals become present in the 

6 see, http://www.eso.org/sci/facilities/paranal/instruments/ 
naco/inst /filters. html for transmission curves of the NACO filters. 



vicinity of the object - not at the object's location itself 
- if we cancel out its peak flux completely. For the rest 
of the analyses and discussion we use AL' ss 9.0 mag 
as default value. The key points of the discussion and 
conclusions remain unchanged for a smaller value of the 
contrast. 

The observed flux densities for HDI00546 based on 
the spectra mentioned above translate into an apparent 
brightness of L' «4.1... 4.3 mag, which fits well to the 
L= 4.02 ± 0.06 mag reported in de Winter et al. (2001). 
Hence we derive an apparent magnitude of L'=13.2 ± 0.4 
mag for the point source component. Here the error is 
the root of sum of squares of 0.3 mag for the uncertainty 
in the contrast and ~0.25 mag for the uncertainties in 
the MIR continuum. Compared to these uncertainties 
the variability of the Pf 7 line flux and the intrinsic error 
in our photometric observations are negligible. 

3.2. Estimating the minimum luminosity 

Assuming that the flux of the point source peaks in 
the L' filter we can estimate its blackbody temperature 
using Wien's law. We can then derive a lower limit on 
the object's luminosity by taking into account its appar- 
ent L' magnitude and its distance. Integrating over all 
frequencies this exercise yields a minimum luminosity of 
L>4-lO- 4 £ . 

3.3. Interaction with the circumstellar disk? 

The VLT/NACO PDI data presented in Quanz et al. 
(2011b) have sufficient spatial resolution and inner work- 
ing angle to probe the disk surface on scales relevant for 
the APP dataset. Those NIR observations revealed sub- 
structures in the disk in the inner few tens of AU. In 
particular the existence of a disk "hole" was suggested 
as both the final polarization intensity images as well as 
the polarization fraction images in H and K s revealed a 
local flux deficit at the same location. 

In Figure 2 we show the large scale disk environment 
revealed by HST/ACS (Ardila et al. 2007) and then, 
zooming in the inner disk regions, the polarization frac- 
tion image of the PDI study. We overplot in red the 
contours of the object detected here. The disk is in- 
clined by ~47°±3° and the position angle of the disk 
major axis is ~138°±4° (Quanz et al. 2011b). If the 
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Fig. 2.— The HD100546 disk on different scales. In the HST/ACS image obtained in the F814W filter (left) the circumsteller disk 
around HD100546 can be traced out to a few hundred AU in scattered light (Ardila et al. 2007). The inner disk regions (~1" in radius) 
are hidden behind the coronagraph or suffer from PSF subtraction residuals. The polarization fraction image (left) obtained at the VLT 
in PDI mode in the H band (Quanz et al. 2011b) probes regions very close to the star, enabling the detection of disk asymmetries not 
accessible with other imaging techniques. The position of the planet candidate is overlaid in the PDI image. North is up and east to the 
left in both images. 



disk surface was smooth and azimuthally symmetric, the 
disk image shown in Figure 2 should be mirror symmet- 
ric with respect to the disk minor axis running with a 
position angle of —48° through the image center (Quanz 
et al. 2011b). However, there are clear asymmetries in 
form of a deficit in polarization fraction in northern di- 
rection, i.e., along the position of the detected object. 
Based on Figure 2 the disk "hole" extends to larger sep- 
arations and appears more like a "wedge" . As discussed 
in Quanz et al. (2011b) the underlying physical reason 
for this feature is not clear at the moment (e.g., drop 
in surface density, disk surface geometry, changing dust 
properties). However, finding an asymmetry at this spe- 
cific location renders plausible a physical link between 
those structures and the source detected here. 



4. DISCUSSION 



4.1. 



The image of an embedded exoplanet? 

Based on the object's position angle, the disk incli- 
nation and the distance to HD100546, the object's de- 
projected separation from the central star is — 68±10 AU, 
i.e., within the large circumstellar disk. Different scenar- 
ios to explain both the L band emission and the observed 
disk structure can be assessed: 

Background source: A background source would be ob- 
served through the HD100546 disk. Based on the disk 
model presented in Mulders et al. (2011) background 
flux in the L band should be attenuated by a factor of 
~6.7-10~ 3 « 5.4 mag at a location of -70 AU 7 . Taking 
this factor into account we used the Besancon galactic 
model (Robin et al. 2003) to estimate the number of ob- 
jects in the apparent magnitude range 7 mag < L < 9 
mag. This yielded —330 objects in a 2 square degree 
patch on the sky centered around HD100546. This num- 

7 This factor does not include that the object is seen through an 
inclined disk which would yield an even higher optical depth. 



ber translates into a probability of having such a physi- 
cally unrelated source in a 1" x 1" field of view around 
the star of p — 1.3 • 10~ 5 . Furthermore, the fact that the 
L band emission appears to be extended argues against 
a background object. 

Disk feature: The observed L' brightness and mini- 
mum luminosity are difficult to explain with disk-internal 
processes alone as the expected temperature in the disk 
mid-plane at the location of the source is only ~50 K 
(Mulders et al. 2011). Furthermore, we are not aware 
of shock-processes that act only locally and might lead 
to the observed luminosity in a disk that appears to be 
not very massive. If it was scattered light that we see, 
one would expect that also in the NIR a maximum in 
scattered light would be seen. Using the PDI images as 
tracer for scattered light we find a local minimum here 
as described above. 

Photospheric emission: If the observed point source 
flux arose solely from the photosphere of a young object 
the COND and DUSTY models suggest masses between 
~15 - 20 Mj U pit or for an age of 5 - 10 Myr (Baraffe 
et al. 2003; Chabrier et al. 2000). Models with lower 
specific entropy in the initial conditions for the formation 
process predict even higher masses (cf. Spiegel & Burrows 
2012). Classical binary formation via core fragmentation 
or formation via disk instability when the disk was still 
massive would be the preferred formation mechanisms 
for an object of this mass. In this case the object formed 
roughly coeval with the star and would have had time 
to significantly alter the structure of the main disk, e.g., 
dynamically clearing a large azimuthal gap, which has 
not been observed. 

Ejected planet: Another massive planet is thought to 
be orbiting in the inner disk gap (e.g., Acke & van den 
Ancker 2006; Tatulli et al. 2011) and we speculate that 
dynamical interactions between multiple planets and the 
disk could have led to an ejection event. The emission 
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we see in the L' images would then be a combination 
of the planet's intrinsic emission plus extra luminosity 
from disk material being heated from the planet moving 
through the disk. Assuming that the planet was initially 
orbiting at 10 AU its orbital period was ^20 years yield- 
ing an orbital velocity of ^3.2 AUyr^ 1 (~15.2 kms -1 ). 
If the ejection velocity is a few times that value it would 
have taken the planet less then 20 years to reach its cur- 
rent location and within less than 100 years it would be 
beyond the extent of the observable disk. Given the age 
of the system, this timescale is extremely small and ob- 
serving the object exactly at the right time is unlikely. 
Adding further complexity to this scenario, the ejection 
needs to occur roughly in the plane of the disk to make 
a link to the observed disk structures. 

Forming planet: In our view the best explanation for 
the observed morphology of both the disk and the emis- 
sion source is the detection of a planet during its forma- 
tion process. The luminosity of the object is not coming 
from an isolated photosphere, but rather the planet is 
still accreting material from the disk. Young, forming 
gas giants with masses between one and a few Jupiter 
masses are expected to have luminosities between 10 — 
10~ 2 £ Q during the first few hundred thousand years after 
gas runaway accretion sets in (Mordasini et al. 2012), in 
agreement with our lower limit. Furthermore, an object 
in this mass range is expected to affect the disk struc- 
ture much less and an azimuthal gap - if it exists - might 
be below our detection limits in the PDI data. A narrow 
gap is hard to see in an inclined disk. This scenario could 
also explain the extended component of the emission de- 
tected here with some disk material being heated in the 
accretion process similar to the case of LkCal5 b 

(see Kraus & Ireland 2012, for possible mechanisms to 
heat the surrounding material during the accretion pro- 
cess). However, we acknowledge that from a theoretical 
perspective the formation of a gas giant planet at this 
location is not readily explainable using first principles. 
For core accretion the timescales to assemble a massive 
rocky core seem to exceed the estimated age of the star, 
and given the observed disk parameters the disk does not 
seem to be gravitationally unstable. 

4.2. Observational tests to distinguish the scenarios 

The different dynamics involved in the scenarios out- 
lined above and multi- wavelength photometry and/or 
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spectra will eventually help us to confirm which hypoth- 
esis is correct. The proper and parallactic motion of 
HD100546 will allow us to rule out a background source 
with second epoch observations obtained as early as mid 
2013. To distinguish between the ejection scenario and 
the "forming planet" scenario, the baseline for follow-up 
observations needs to be a few years. While the object 
is expected to orbit its star with a period of ^360 years 
in the latter case, it should move away quickly in radial 
direction if it were ejected. Also, high spatial resolution 
ALMA observations will help to search for an azimuthal 
gap in the surface mass density (gas and/or dust) at the 
planet's location. Spatially resolved information about 
the existence (or non-existence) of an azimuthal disk gap 
may allow us to derive a dynamical mass estimate for the 
planet (Lin & Papaloizou 1986; Bryden et al. 1999). 



5. CONCLUSION 

We have presented observational evidence that a gas 
giant planet could be forming in the circumstellar disk 
around the Herbig Ae/Be star HD100546 at a separa- 
tion of ^68 AU. This scenario, among others that we 
discussed, seems best capable of explaining most of the 
available data. However, some aspects remain qualita- 
tive and follow-up observations are required to validate 
our proposed interpretation. Together with LkCal5 b 
(Kraus & Ireland 2012), the object presented here is cur- 
rently the best candidate for a forming young gas giant 
planet. Particularly interesting is that HD100546 "b" 
would be the first protoplanet that is still embedded in 
a circumstellar disk and that it forms at large orbital 
separations. If confirmed, HD100546 would be a unique 
laboratory to study planet formation and the interaction 
between forming planets and the disk directly. 
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